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Abstract 

We study mechanism of formation of synchronized clusters in coupled maps 
on networks with various connection architecture. The nodes in a cluster 
are self synchronized or driven synchronized, based on the coupling strength 
and underlying network structures. Smaller coupling strength region shows 
driven clusters independent of the network rewiring strategies, where as larger 
coupling strength region shows the transition from the self-organized cluster 
to the driven cluster as network connections are rewired to the bi-partite 
type. Lyapunov function analysis is performed to understand the dynamical 
origin of cluster formation. The results provide insights into the relationship 
between the topological clusters which are based on the direct connections 
between the nodes, and the dynamical clusters which are based on the func- 
tional behavior of these nodes. 
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1. Introduction 

Networks of dynamical units have been used to model everything from 
earthquakes to ecosystems, neurons to financial market 
Studying the impact of complex network topology 0, 0] on the dynamical 
processes is of fundamental importance for understanding the functioning 
of these real world systems. Synchronization is one of the most fascinating 
phenomenon exhibited by such nonlinear dynamical units interacting with 
each other [§[ 0]. This phenomenon has been intensively investigated as a 
universal concept in many disciplines, particularly in neuro-sciences [10!] and 
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communication 



111 ] . Recently synchronization has been shown to play a 



crucial role in the species extinction [12] . Keeping such applications in mind, 



a great deal of work is devoted to build up the strategies for the network 
connection architecture which could synchronize the dynamical units on them 



131 ] . Most of these works have considered global synchronized states spanning 
all nodes in a network. 

In the past years focus of coupled dynamics research got shifted to the 
(phase) synchronized clusters [l4j or dynamical clusters instead of a global 
synchronized state. Dynamical clusters are defined based on some kind of 
coherent time evolution of some units which is different than the evolution 
of other units. Dynamical clusters have been the part of coupled dynam- 
ics research since the seminal work of Kuramoto which mainly focused on 



the transition to the global synchronized state [15j. Before the transition, 
(loT ] reports the phase- synchronized clusters, where clusters are distinguished 
based on the frequencies of dynamical units. Kaneko's coupled maps model 



also talks about the glassy phase where nodes are partially synchronized [16 



Kuramoto's coupled oscillator and Kaneko's coupled maps models assumed 
either nearest neighbor connections or all to all connections between the 
nodes. Later on with the spurt of complex network research, these types of 
coupled dynamics models had been investigated with units having local and 
nonlocal connections, and also show partial synchronized state with dynam- 



ical clusters [17(. These works mainly concentrated on the finding of partial 
synchronized state, and did not pay much attention to the coupling con- 
figurations. However, networks modeling real world, have some structures, 
hence it was very natural to imagine a relation between the connection ar- 
chitecture and the dynamical clusters. It could be rephrased as follows: do 
connections in the dynamical clusters have some specific features ? Man- 
rubia and Mikhalov in [18] found that the connections in dynamical clusters 
do have some particular configuration, i.e. elements from a cluster gener- 
ally having more connections inside the cluster than with the elements from 
other dynamical clusters. Though this was a very significant observation 
and more importantly presented the concept of relation between dynami- 
cal cluster and the underlying network, their observation was based on the 
dynamical units coupled with some particular coupling strength and having 
random coupling architecture. The detailed investigations of coupled dynam- 



ics on various networks, done by Jalan and Amritkar [19[, revealed that there 



exists two different mechanisms of the formation of the dynamical clusters; 
(a) self-organized clusters where couplings are intra-cluster type, (b) driven- 



2 



synchronized where coupling are inter-cluster type. The main result of the 
Ref . [19| is that the coupled dynamical units may get synchronized because of 
inter-cluster or intra- cluster couplings. Whereas self-organized mechanism 
has more commonly been thought as a reason of synchronization from semi- 
nal work by Kaneko 16] to recent work by several scientists 



18, 20, 21 , 1 



identified a new mechanism of cluster formation that is driven synchroniza- 
tion and did extensive studies of the two mechanisms of cluster formation for 
various networks 19|, |22|, |23j . Recently self and driven-synchronization has 
been investigated in the relevance of brain cortical networks [24| . 

Real world networks have community or module structure [251.1261]. where 
modules are the division of network nodes into various groups within which 
the network connections are dense, but between which they are sparser. The 
modularity concepts assumes that system functionality can be partitioned 
into a collection of modules and each module is a discrete entity of several 
components performing an identifiable task, separable from the functions of 
other modules 26] . In the light of this finding the concept of relation between 



dynamical clusters and the network architecture [18|, [19| can be rephrased; 
do dynamical clusters reveal organization or module structure of the under- 
lying network, and what functional clusters would be preferred for a given 



ie_p 

network architecture ? Investigation in this direction 27|, |28| showed that the 
dynamical units reveal the hierarchical organization of the synchronization 
behavior with respect to the collective dynamics of the network. Furthermore 
Zhou et. al. used synchronized clusters to unveil the functional connectivity 
in the complex brain networks j29|. These works exploited the synchronous 
behavior of dynamical units to detect the hierarchical organization, on the 
other hand M. Timme concentrated on the disordered state to understand 
the relation between dynamical behavior and the network topology 30] . Re- 



cently phenomenon of driven synchronization is observed for three stochastic 



oscillators which are coupled with delay [31 



In this paper we highlight that there is a great difference in the process 
leading to the synchronization and synchronized cluster state, depending 
upon the particular coupling configuration. We show that these two ways of 
the construction of the dynamical clusters depend highly on the underlying 
network architecture. Self-organized clusters are formed for the networks 
whose network architecture is such that it can clearly be divided into dif- 
ferent sub-communities, whereas for the bi-partite networks driven clusters 
are formed. Mixed clusters, i.e. some are self-organized and some are driven 
organized, are formed for the network architecture which lie in the middle 
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of these two extreme ideal network structure. We use Lyapunov function to 
understand the origin of these two types of clusters. 

The paper is organized as follows. After the introductory section describ- 
ing the current understanding and importance of dynamical cluster research, 
we introduce general model of coupled dynamics in section II. Section III 
consists the definition of dynamical clusters and measures quantifying the 
behavior of the clusters. We explain network rewiring strategy in section IV 
and also present results for the dynamical behavior at various step of the 
network rewiring. Section V considers Lyapunov function analysis. Section 
VI concludes the paper with the possible future directions. 



2. Coupled dynamics on Networks 



Consider a network of N nodes and N c connections. Let each node of the 
network be assigned a dynamical variable Xi, i = 1, 2, . . . N. Evolution of the 
dynamical variable is written as 



16 



Xi(t 



f(xi(t)) + ^ C ^9Mt))-9^M)} 



(1) 



where Xi(t) is the dynamical variable of the ith node at the tth time step, 
and C is the adjacency matrix with elements taking values one and zero 
depending upon whether i and j are connected or not, s G [0,1]. Matrix 
C is symmetric matrix representing undirected network, and ki = J? . CV,- 
is the degree of node i. Function f(x) defines the local nonlinear map and 
function g(x) defines the nature of coupling between nodes. In this paper we 
present results for the local dynamics given by logistic map f(x) = fix(l — x) 
in the chaotic regime and g(x) = f(x). This coupled maps model is studied 



extensively to get insight into the behavior of real systems [la . |32j. We 
take the value of // = 4, for which individual un-coupled unit shows chaotic 
behavior with the value of Lyapunov exponent being In 2. As an effect of 
coupling the coupled dynamics [1] shows various different kinds of coherent 
behavior, such as synchronization [l6j], phase-synchronization 14, 19, 22, 
33] and other large scale macroscopic coherence [35| depending upon the 
coupling architecture and the coupling strength. In this paper we consider 
synchronization Xi(t) = Xj(t) for all i, j and all initial conditions, and phase- 



synchronization as considered in [3J 
of times the variables Xi(t) and Xj(t) 



19 



t 



Let rii and rij denote the number 
= 1, 2, r for the nodes i and j 



show local minima during the time interval r . Let rij,- denote the number 
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of times these local minima match with each other. Phase distance between 
nodes i and j is then given by dij = 1 — 27iy/(n* + rij). Node i and j are 
phase synchronized when = 0. In a dynamical cluster, all the pairs of 
nodes are (phase) synchronized. 



3. Synchronized clusters 

Various types of dynamical clusters based on the couplings between the 
nodes are defined as: (a) Self-organized cluster: If nodes of a cluster are 
synchronized because of intra-cluster couplings, the cluster is referred as self- 
organized cluster, and (b) driven cluster: If nodes of a cluster are synchro- 
nized because of inter-cluster connections, the cluster is referred as driven 
cluster. Quantitative measure of above two types of clusters are given by the 



following two quantities [19l . |22 



f intra 
f inter 



N, 



intra 



N, 



inter 



where Ni ntra and Ni nter are the number of intra- and inter cluster connections, 
respectively, f inter and f intra would take values between zero and one. For 
the ideal driven or ideal self-organized case one of the measure would take 
value one. According to this definition j : inter and f intra both would be zero 
when none of the nodes form cluster. The intermediate values of these two 
quantities represent the situation where some of the nodes are of driven type 
and some of them are of self-organized type, or it may correspond to the 
case when some clusters have mixed connections, i.e. some of connections 
are intra-cluster type whereas some are inter-cluster type, f inter + f inter = 1 
represents the situation when all the nodes form one or many synchronized 
clusters and < f intra + f inter < 1 corresponds to the case when some nodes 
do not belong to any cluster and evolve independently. In the dominant self- 
organized clusters number of intra-cluster connections are much greater than 
the inter-cluster connections, while in the dominant driven clusters reverse 
is true. Mixed clusters corresponds to the case when number of inter and 
intra-cluster connections are roughly equal. Note that while counting f inter 
we do not consider isolated nodes, only the nodes which are the part of any 
synchronized cluster are considered. 
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Newman's definition of modularity of a network, Q, is also based on the 
ratio of intra and inter connections (Q = l/2iV c [£\e# — Y^ijk e ij e ki))i where 
Cij is the fraction of edges in the network that connect nodes in group i to 
those in groups j. Modularity of a network is one if it can be divided into 
subgraphs such that nodes in each subgraphs are connected within it, with 
a single connection to any of the other subgraph. At least one connection is 
required in order to keep the network connected. We would call these clus- 
ters, structural clusters (or community) to distinguish it from the dynamical 
clusters. In self-organized state, nodes in same community get synchronized 
dynamically also, while for the driven state, nodes which are synchronized 
have connections only to the other dynamical clusters. 

In the first case, it is more obvious, that the nodes which are connected, 



tend to synchronize with one another. This concept was exploited in |29| to 



find out functional hierarchy in cat's cortex network. However, as shown in 



the Refs. 19|, |22|, |23|], dynamical clusters could be formed without having 



even a single connection within the cluster (ideal driven state). These two 
types of clusters, self-organized and driven, could be observed irrespective of 
the dynamical state of the EqJTJ i.e. whole dynamics may lie on chaotic, pe- 
riodic, quasiperiodic attractor depending upon the coupling strength and the 
connection architecture. Different possible states of this coupled dynamics 



are discussed in details in 22]. This paper focuses on the relation between 
different types of clusters with the structure of the coupling matrix C. We 
particularly track the nature of dynamical cluster with the rewiring of net- 
work and show that depending upon the modularity of network, the type of 
dynamical clusters changes. If network is completely modular then in general 
self-organized clusters are preferred, whereas for bi-partite kind of network 
driven clusters are preferred. 

In the following we present results for the coupled maps on networks, as 
described by the Eq. (pQ), evolved with random initial conditions. First the 
initial network and the rewiring procedure are explained and then various re- 
sults demonstrating the relation between the dynamical clusters and network 
structure are demonstrated. 



4. Evolution of synchronized clusters with Network rewiring 

We start with a network of size N having two complete subgraphs, G\ and 
G2 with N\ and N2 nodes respectively. Complete subgraph means, all the 
nodes in G\ (G2) are connected with the all other nodes in G\ (G2) (Fig. Ufa)). 
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Figure 1: Pictorial representation of networks at different rewiring step, (a) Initial network 
with two complete subgraphs G\ and G2 (b) rewired network at n = 1, solid lines are the 
original connections and dotted lines are the new connections which a node from the group 
Gi makes with the nodes of G2, (c) final bi-partite graph. 




1 

- 


1 




\ fj (b) _ 


_ a 









1 

- 


1 1 1 1 








/*• — v 








1,1,1, 






ft 1 1 






1 * 








(f) 



Figure 2: (Color online) Fraction of inter (•) and intra (o) connections as a function of the 
coupling strength e. After an initial transient (about 2000 iterates) phase synchronized 
clusters are studied for r = 100. The logistic map parameter is = 4, and Ni — N2 = 50. 
(a) for the original network at rewiring step n = 0, (b) n = 5, (c) n — 10, (d) n = 20, 
(e)n = 30, and (f)n = 50. All figures are plotted for the average over 20 sets of random 
initial conditions for the coupled dynamics. 
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Figure 3: (Color online) The fraction of number of clusters M c /„ s (o) and the fraction of the 
number of nodes forming clusters N c i ua (») as a function of the coupling strength e. The 
network structure are same as for the Fig. (2) and quantities are plotted for the average 
over 20 realization of random sets of the initial conditions. 



There exists one connection between G\ and G 2 , which is necessary to keep 
the whole graph connected. Now at each rewiring step (n) disconnect one 
node from subgraph G\ (G7 2 ) and connect it to all the nodes in the subgraph 
G2 {Gi) (Fig.[D^b)). With this rewiring scheme, average degree of the network 
remains of the same order (N). After rewiring step Nx, all nodes in G\ are 
connected with all the nodes of G 2 , and for the case Nx = N 2 , the graph is 
complete bi-partite (Fig. [U(c)). At each step of the rewiring we evolve Eq. (pQ) 
with random initial conditions and study the nature of the dynamical clusters 
after some initial transient. For small coupling strengths region (e < 0.5), 
number of nodes forming synchronized clusters is very small and hence in 
this region we consider phase-synchronized clusters. For the larger coupling 
strengths, in general, nodes form exact synchronized clusters. In the following 
we consider dynamical clusters based on the phase-synchronization. Fig. (J2]) 
shows behavior of f intra and J \ nter as a function of coupling strength e for the 
various steps of the network rewiring according to the above strategy. 

Fig. ([2]) shows that after initial turbulent regime where none of the node 
synchronizes , nodes form synchronized clusters. For the weak coupling (e ~ 
0.2) the coupled dynamics governed by Eq.[T]shows partial synchronized state 
corresponding to the many dynamical clusters with finter and f intra both 
being nonzero. Subfigure (a) is plotted for the initial network (Fig. 11(a)), 
subfigure (b) is plotted for the rewiring stage n — 5, which means that 5 
nodes from each group are rewired. These nodes break the connections with 
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their communities and make new connections with the nodes of the other 
communities. 

Fig. (j3J) plots the fraction of number of clusters (Ma us ) and the fraction 
of number of nodes forming clusters N c i us . The first measure M c \ us counts 
an isolated node as a separate cluster, and the second measure N c i us counts 
only those nodes which are the part of a cluster and are not isolated. It is 
clear from the subfigures that for each rewiring step almost all the nodes 
form cluster after some coupling strength e > 0.2 value. Note that for the 
lower coupling strength region, in general, nodes form phase- synchronized 
cluster, as defined in the previous section, whereas for the stronger coupling 
strength region e > 0.5 they form exact synchronized clusters. For our 
investigations it does not matter whether nodes are phase-synchronized or 
exact synchronized, as long as large number of nodes form the clusters we 
get relevant information regarding dynamical cluster and network structure. 
It is seen from all the subfigures of Fig. (jSJ), that the nature of synchronized 
clusters is very much similar in the weaker coupling strength regime (e ~ 
0.2) for all rewiring steps, with (a) and (f) being the extreme cases of the 
network structure. Hence, for our study the weaker coupling strength region 
is not very interesting as in this region one does not seem to get relation 
between dynamical cluster and network structure. In order to get insight 
about the relation between dynamical clusters and the network structure 
we concentrate on the coupling strength region e > 0.5. Coupled dynamics 
Eq. ([I]) gives interesting dynamical clusters in this coupling strength region, 
where connections between the nodes of dynamical clusters differ as structure 
of the network is varied. In the following, the behavior of coupled dynamics 
and nature of synchronized clusters for the various coupling strengths in this 
region are discussed. 

Fig. 01]) plots the synchronized clusters for different networks generated 
with the above rewiring strategy. The figures are plotted for the coupling 
strength e = 1. For this coupling strength the nodes form exact synchro- 
nized cluster rather than the phase-synchronized cluster which is the case for 
weaker coupling strengths. Fig. (@J)(a) is for the rewiring step n = 0, when 
nodes in the network form two complete subgraphs. The dynamical behavior 
of these nodes show synchronization with two clusters, and these clusters are 
of self-organized type; i.e. nodes belonging to the structural clusters form 
dynamical clusters. Fig. (@J(b) is node-node plot for the network at step 
n = 20, i.e. connections of 20 nodes from each group are rewired such that 
the 20 nodes from the group G\ (G2) get connected with the 30 nodes of 
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Node Node 

Figure 4: (Color online) Node vs node diagrams illustrating the behavior of nodes at 
different step of rewiring. After an initial transient synchronized clusters are studied. 
The solid circles (•) show that the two corresponding nodes are coupled and the open 
circles (o) show that the corresponding nodes are (phase) synchronized. In each case the 
node numbers are reorganized so that nodes belonging to the same cluster are numbered 
consecutively. All plots are for coupling strength e = 1. (a) for the network having two 
complete subgraphs with Ni — N2 = 50 (b) network at the rewiring step n — 20 (c) for 
n = 30 and (d) for n — 50 which leads to the complete bipartite network. 

the group G2 At this stage mixed dynamical clusters are formed. Two 

small groups of nodes which are rewired to the different groups loose the 
synchronization with the nodes in their respective group and are synchro- 
nized independently forming two separate clusters. These two clusters, first 
and fourth cluster from the bottom, are therefore of driven type. Rest of the 
nodes, which remain fully connected inside their respective groups, remain 
self-organized type. Fig. (jSJ) shows dynamical clusters at the various stage 
of rewiring for the coupling strength value e = 0.9. Subfigure (a) shows the 
two self-organized clusters for the initial network (Fig. [p. Subfigures (b), 
(c) and (d) are plotted for the n = 20, 30 and 50 rewiring steps respectively. 
For n = and n = 50, self-organized and driven clusters respectively are 
the stable configurations, thus coupled dynamics lead to these for both the 
coupling strengths. For intermediate rewiring states dynamical clusters vary 
with the value of coupling strengths. Node numbers are reorganized so that 
nodes belonging to the same cluster are numbered consecutively. Fig. ([6]) 
plots largest Lyapunov exponent as a function of the rewiring step n, n = 
corresponding to the original network (Fig. QJa)). At each rewiring step one 
node from each group is rewired. After n = Ni = N 2 = 50, the network 
becomes like Fig. CQ)(c). After fixing the coupling strength, Lyapunov expo- 
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100 



Node Node 

Figure 5: (Color online) Node vs node diagram illustrating the behavior of nodes (same 
as Fig. (4)), for the coupling strength e = 0.9. Network structure remains same as of the 
Fig. (4). Node numbers are reorganized so that nodes belonging to the same dynamical 
cluster are numbered consecutively. 




Figure 6: (Color online) Largest Lyapunov exponent A as a function of the rewiring step 
n. Figures are plotted for two different coupling strengths e = 0.9 and e = 1. Lyapunov 
exponents are calculated for the average over 20 realization of the random set of initial 
conditions. 
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nent A is: calculated for the coupled dynamics (Eq. [Q at each step of the 
rewiring. Results for e = 0.9 and e = 1 show that coupled dynamics may 
lie on the periodic (A < 0), quasiperiodic (A ~ 0) or on the chaotic (A > 0) 
attractor, but the nature of the coupled dynamics governed by the Eq. [I] 
remains as shown by the Fig. (TJ]) andO depending upon n. 



5. Lyapunov function analysis 



We define Lyapunov function for any pair of nodes i and j as [36|, l37| 

Vy(t) = (Xt(t) - Xj (t))* ( 2 ) 

Clearly, V^-(t) > and the equality holds only when the nodes i and j are 
exactly synchronized. For the asymptotic global stability of the synchronized 
state in a region, Lyapunov function should satisfy the following condition 
in that region: 

vit + i) 

V(t) 

For the global synchronous state (xi(t) = Xj(t), Vi, j), we can write the con- 
dition for two synchronized clusters in the self-organized state (i.e. %i x if) = 
x i2 (t); Via , «2 ^ Gi an d Xj^t) = x j2 (t); Vji, j 2 e G 2 ), as; 



V hi2 {t+l) = [(l-e)[f(x k (t))-f(x i2 (t))} 

[fixtM-fixM)} 



iVi - 1 



(3) 



Using Eq. [3] and the above equation for coupled dynamics (Eq. [T]), we get the 
coupling strength region for which the synchronized clusters state is stable; 

iVl - 1, 1, 

— (!--)<«<! (4) 

With the rewiring, say at rewiring step n = u, we get v nodes from each 
group rewired such that now there are four different types of nodes. One can 
quickly see by using Lyapunov function test that the following synchronized 
clusters state would be stable ; 

1) Cluster one(Si) with the nodes which remain unwired in the group G\, 
synchronized dynamics of nodes in this cluster is Xi(t) = Xi(t). 

2) cluster two(S2) with nodes from group G\ which get disconnected with 
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Figure 7: Dynamical clusters at a rewiring state n t . Si , S2 , £3 and S4 are different dy- 
namical clusters as described in the text. G\ and Gi are the dynamical clusters for the 
initial network n = 0. 



the other nodes in G± and get connected with all the nodes in G 2 , %i = X^it). 

(3) cluster three(Sz) with nodes which remain in group G2, Xi = Xs(t). 

(4) cluster four(Si), rewired nodes in G 2 , Xi(t) = X 4 (t). 

Lyapunov function for the nodes in the cluster one (and for the cluster three) 
remains same as EqJHJ because coupling of the nodes in this group to the 
groups three and four (one and two) cancel out. Lyapunov function for the 
nodes in group 5*2 (S4) can be written as; 

V^(t + l) = [(l-e)(/(x i (t))-/(x i (t)))] a 

This equation is very simple because the coupling terms for the nodes % and 
j are same and hence they cancel out. Condition for this state to be stable 
is given as; 

Vv(t) ~ 1 } 

+0{{x l {t)-x ] {t)) 2 )] 2 <1 

For f(x) = fix(l — x) and using < Xi(t) + Xj(t) < 2 and condition for 
achieving synchronized state, we get the range of coupling strength as ((/i — 
1)///. < e < 1) for which the synchronized cluster S2 (as shown in the Fig. (J7J) 
is stable. Similarly we can write the condition for the nodes in the group 
G2 which are rewired and make synchronized cluster as S4, and nodes which 
make synchronized cluster as S3. 

6. Conclusions and Discussions 

In order to explore the relations between dynamical clusters and net- 
work clusters we study coupled maps on various networks generated with a 



f\xm + Xi{t) 2 Xj{t) f'^At)) 
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simple rewiring strategy. Starting with a network having two complete sub- 
graphs, nodes are rewired at each step such that after some rewiring steps we 
get a complete bi-partite network. Rewiring strategy is adopted such that 
average degree of the networks at each rewiring step remains of the order 
of N. Smaller coupling strength values show phase synchronized clusters of 
dominant driven type, whereas larger coupling strength region show different 
mechanisms of synchronized clusters formation depending upon underlying 
network architecture, and hence provide insight into role of network archi- 
tecture in the coherent behavior of the associated dynamical units. Coupled 
dynamics form self-organized clusters for the network having two complete 
subgraphs, and form driven clusters for the bi-partite case. At intermediate 
steps, the nodes receiving similar input form dynamical clusters, and these 
clusters could be self-organized type, driven type or mixed type depending 
upon which connection environment they belong to. Lyapunov function anal- 
ysis shows that for the driven synchronization if any two nodes have similar 
coupling architecture, the difference of the dynamical variables for these two 
nodes cancel out. Whereas for the self-organized synchronization, the cou- 
pling term corresponding to the direct coupling between the nodes do not 
cancel out, and other couplings which are common to both the nodes cancel 
out. The 



Ref. [23j shows that there are two mechanisms of cluster formation in 
networks, and using the global stability analysis it provides arguments be- 
hind the mechanisms by taking globally coupled and complete bipartite as 
the extreme cases. Here we show that there is a gradual transition from the 
self-organized to the driven behavior, as the underlying network is rewired 
from the two clusters to the bipartite network. The small coupling val- 
ues do not show any impact of the structural changes on the mechanism of 
the synchronization, whereas large coupling values show significant signature 
of the underlying structure on the synchronized clusters. The dense net- 
works (N c ~ N 2 ) considered here yield all the nodes forming clusters at each 
rewiring stage, and clusters are of different types depending upon the under- 
lying network structure at that particular rewiring step. Through extensive 
numerical simulations, we show that the nodes, getting similar coupling en- 
vironment, are synchronized irrespective of whether they are connected or 
not. At the intermediate rewiring steps and for the higher coupling values, 
few clusters are mixed type, and few clusters are of the ideal driven or of 
the ideal self-organized type. All the nodes in each cluster receive similar 
coupling environment. Using the Lyapunov function analysis, we get a clear 
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picture of the synchronization for each cluster individually. The stability 
condition for the synchronization of the nodes in clusters S^ and S4 leads to 
the reasoning valid for the ideal driven cluster. Also, stability range written 
for the dynamical clusters Si and S3 matches with that of the ideal self- 
organized cluster because of the very simplistic underlying network model 
we are using, where any terms outside these cluster cancel out leading to the 
condition shown by the globally coupled network in the synchronized state. 

Though the network architecture and rewiring strategies considered here 
are very simple, whereas real world networks have complicated random struc- 



ture or have complicated rewiring or connection evolution strategies 38], the 
results presented here shows direct impact of network connection architec- 
ture on the evolution of dynamical units. The main aim of the analysis here 
is to emphasize that the dynamical clusters do have information about the 
structure, but the formation of driven clusters reveal that the nodes in the 
dynamical clusters are not always the nodes which are connected, rather the 
nodes which have similar environment show coherent behavior and form a dy- 
namical cluster. For a given network, depending upon the coupling strengths 
there could be various possible states which form stable configurations of the 
dynamical clusters. Results presented in the paper indicate that synchronous 
behavior of the dynamical units interacting with each other might be receiv- 



m a 



ing the similar inputs, for instance synchronous firing of neurons [10 
cortex may be because of the couplings among the neurons, or may be be- 
cause of getting similar synaptic inputs from the neurons of different cortical 



areas. In the continuation of previous work |19|, [22|, [24| , the results presented 
here strengthens understanding about the possible relations between the dy- 
namical behavior and the underlying connections architecture. Future inves- 
tigations would involve more general network rewiring or evolution scheme, 



as well as effect of delays |39[ on the phenomenon of driven synchronization. 
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